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The localization of a quantum state is numerically exhibited in a nonuni- 

> 

r**"- ■ tary Newtonian model for gravity. It is shown that an unlocalized state of a 

ball of mass just above the expected threshold of 10 11 proton masses evolves 

O 

(T) . into a mixed state with vanishing coherences above some localization lengths. 



03.65.-w, 03.65.Bz 



3 ' I. INTRODUCTION. 

cr 

The apparently inescapable conclusion of the so called "information loss paradox" was 
the assumption that a full theory of quantum gravity would be non-unitary, with pure 
states evolving in principle into mixed ones [1]. At the basis of such assumption was the 
observation that the area of the event horizon of black holes increases monotonically in time 
while information is irretrievably lost when matter falls into black holes. This in turn led 
Bekenstein to propose that a generalized second law holds, where the entropy of a black 
hole is proportional to the area of its event horizon [2]. More generally the laws of black 
hole dynamics exhibit a remarkable resemblance with the laws of thermodynamics, if energy 



*Unita I.N.F.M., I.N.F.N. Salerno 



1 



is replaced by the black hole mass and temperature is assumed to be proportional to the 
surface gravity. This idea was dramatically confirmed by Hawking's discovery that black 
hole evaporation has a thermal character with a temperature proportional to the surface 
gravity [3]. 

It was argued that a non-unitary evolution should be effective even in ordinary labora- 
tory physics, since a continuous creation and evaporation of submicroscopic black holes is 
expected within such a theory, even though it seemed to give rise to unacceptably large vio- 
lations of causality or energy-momentum conservation [4,5]. A more recent analysis pointed 
out that, in general, violations could be undetectable at the scales of laboratory physics, and 
in particular that there need not be any (easily detectable) conflict between non-unit arity, 
causality and energy conservation if the model is non-Markovian, as suggested by black hole 
dynamics [6]. This analysis, though inspired by black hole dynamics, was performed for 
generic non-unitary models without any specific reference to gravity and intended only to 
mimic a possible effective model for gravity at low energy. 

It is remarkable that indications for a non-unitary evolution induced by gravity also 
emerge in the very context of non-relativistic QM, in which many tentative ad hoc modifica- 
tions of the Schroedinger equation ("reduction models") have been proposed to account for 
both wave-function collapse and the emergence of classicality [7-9]. In fact in these models, 
in order to minimize energy pumping one is induced to couple the phenomenological stochas- 
tic external field mainly to the centre of mass, which strongly suggests its gravitational origin 
[10,11]. 

In view of a possible convergence of these two research lines (information loss paradox 
and reduction models), in some recent papers one of the authors proposed and analyzed a 
non-Markovian non-unitary model for Newtonian gravity [12,13]. It has been shown that this 
model can be derived from fourth-order gravity in which, while classical instability is cured 
at expense of unitarity, the nice one-loop ultraviolet behavior of HD gravity is preserved 
[14,15]. 

While for the non-Markovian models considered by Unruh and Wald the basic idea is to 
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have the given system interacting with a "hidden system" with "no energy of its own and 
therefore... not... available as either a net source or a sink of energy" [6], in the present 
model energy conservation is granted by the "hidden system" being a copy of the physical 
system, coupled to it only by gravity, and constrained to be in its same state and then to 
have its same energy. The unitary dynamics and the states referred to the doubled operator 
algebra are what we call respectively meta-dynamics and meta-states, while, by tracing out 
the hidden degrees of freedom, we get the non-unitary dynamics of the physical states. Pure 
physical states correspond then to meta-states without entanglement between physical and 
hidden degrees of freedom. 

In particular it has been shown that, while reproducing the classical aspects of the New- 
tonian interaction, the model gives rise to a threshold, which for ordinary condensed matter 
densities corresponds to ~ 10 11 proton masses, above which self-localized center of mass 
wave functions exist, in analogy to the Schroedinger- Newton model [16-19]. However, while 
the latter produces rather unphysical stationary localized states without spreading, accord- 
ing to the present model an initial localized pure state evolves in time into an unlocalized 
ensemble of localized states. That is consistent with the expectation that (self-) gravity may 
produce a growing entropy in a genuinely isolated system [20,15], as suggested by black hole 
dynamics. 

It was also shown that the meta-dynamics can be reformulated in terms of a functional 
integral over two auxiliary scalar fields. If the resulting expression is applied to a physical 
state, obtained by tracing out the hidden degrees of freedom, one gets an expression which 
can be taken as an independent definition of the model, free from any reference to hidden 
degrees of freedom [13,15]. This expression was used both to show that the N — > oo limit 
of the iV-hidden replicas generalization of the model reproduces the Schroedinger-Newton 
model and to analyze wave function reduction [21]. To be more specific, the superposition 
of a large number of widely spaced localized states, corresponding to the ground meta- 
state of the relative motion of a physical body of mass M and its hidden partner, was 
considered. It was shown that, within a random phase approximation and omitting the slow 
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spreading of the localized states, it evolves into an ensemble of the mentioned localized states 
in a typical time r g = HG -1 M~ 5 / 3 q -1 / 3 ~ 10 20 (M/m p )~ 5 ^ 3 s (for ordinary matter density 
q ~ 10 24 m p / cm 3 ), where G and m p respectively denote the gravitational constant and the 
proton mass. 

The key hypothesis behind any model of emergent classicality is that for some reason 
(the surrounding environment or a fundamental noisy source) unlocalized macroscopic states 
get localized, i.e. that coherences for space points farther than some localization lengths 
vanish [22]. In such a way embarrassing quantum superpositions of distinct position states 
of 'macroscopic' bodies are avoided, thus recovering an essential element of reality of our 
classical realm of predictability: the elementary fact that bodies are observed to occupy quite 
definite positions in space. The aim of the present paper is to check that feature, which makes 
the model a viable localization model, without using very peculiar initial conditions or any 
approximation, and independently of the heuristic functional formulation. In order to make 
the problem numerically tractable we take a rotational invariant initial state of an isolated 
ball of ordinary matter density and a total mass just above the mass threshold, namely in 
the least favorable conditions to exhibit dynamical localization. The ensuing phenomenology 
may be also relevant in view of possible future experimental checks of fundamental (not 
environment induced) decoherence. In fact in the last few years much efforts have been 
done in order to get a quantitative control of the environmental decoherence [23] and to 
construct cat states involving growing masses. In this respect an improvement in the Bose- 
Einstein Condensation (BEC) technology may be useful to test gravity-induced decoherence, 
considering the unprecedented scale of controlled quantum coherence achieved there. 

II. THE MODEL. 

We give here a concise definition of the model. Let H[ifj\ifj] be the non-relativistic 
Hamiltonian of a finite number of particle species, like electrons, nuclei, ions, atoms and/or 
molecules including also the halved gravitational interaction, where tp^,ip denote the whole 
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set tpj(x),ipj(x) of creation-annihilation operators, i.e. one couple per particle species and 
spin component. H[ip\ip] includes the usual electromagnetic interactions accounted for in 
atomic, molecular and condensed-matter physics. To incorporate that part of gravitational 
interactions responsible for non-unitarity, one has to introduce complementary creation- 
annihilation operators ipj(x),ipj(x) and the overall (meta-)Hamiltonian 

1 j,k J |X 2/1 

acting on the product ® of the Fock spaces of the ip and ip operators, where rrii is the 
mass of the i-th particle species and G is the gravitational constant. The ip operators obey 
the same statistics as the corresponding operators ip, while [ip,ip]- = [ip,ip^]- = 0. 

The meta-particle state space S is the subspace of F^ <g> F^ including the meta-states 
obtained from the vacuum ||0)) = |0)^, <8> |0)^ by applying operators built in terms of the 
products i/)j(x)ipj(y) and symmetrical with respect to the interchange tp^ <-> ip^, which, 
then, have the same number of ip (physical) and ip (hidden) meta-particles of each species. 
As for the observable algebra, since constrained meta-states cannot distinguish between 
physical and hidden operators, it is identified with the physical operator algebra. In view of 
this, expectation values can be evaluated by preliminarily tracing out the ip operators. In 
particular, for instance, the most general meta-state corresponding to one particle states is 
represented by 

||/)) = JdxJ dyf(x,y)^(x)^(y) |0> , f(x,y) = f(y,x). (2) 
This is a consistent definition since Hq generates a group of (unitary) endomorphisms of S. 



III. DYNAMICAL LOCALIZATION: NUMERICAL RESULTS. 



Consider now a uniform matter ball of mass M and radius R. Within the model the 
Schroedinger equation for the meta-state wave function E(X,Y,t) is given by 



OS 
dt 



2M 



{V 2 X + V 2 Y ) + V{\X-Y\ 



(3) 



where X and Y respectively denote the position of the center of mass of the physical body 
and of its hidden partner, while V is the (halved) gravitational mutual potential energy of 
the two interpenetrating meta-bodies, which, as can be shown by an elementary calculation, 
is 

yn _ 1 rM 2 (9(2R-r)(80R 3 r 2 -30R 2 r 3 + r 5 - 192i? 5 ) 9(r - 2R)\ 
V{r)--GM ^ - J (4) 

where 9 denotes the Heaviside function. Observe that our final result can be immediately 
reread as the solution for the whole set of parameters obtained by the scaling: 

t^Xt, M — > A~ 1/5 M, X -> A 3/5 X, R -> X 3/5 R 

E (X, Y; t) -> S (A~ 3/5 X, A~ 3/5 F ; t) (5) 

where A is a real positive dimensionless parameter. This is consistent with the expression for 
T g , the latter giving r g — > \r g . Besides, for consistency, we note that the mass cannot cross 
the threshold, as the latter scales with the same power law of the mass itself M t — > A _1 / 5 M t . 

If we separate Eq. (3) into the equation for the relative motion and that for the center 
of mass, then for ordinary matter density g ~ 10 24 m p / cm 3 and M above the threshold 
~ 10 n m p , the former admits bound meta-states of width Ac ~ {m p jM) x l 2 cm, corresponding 
to small oscillations around the minimum of the gravitational potential. In particular an 
untangled localized meta-state, corresponding to a physical pure state is: 

m TOT = ^oPO^oOH = Mi* + Y] /2)* ([X - Y] /2), (6) 

where 

" Y] j2) = (A^)- 3/4 exp[- \X-Y\ 2 /(2A 2 G )}; A G = (8h 2 R 3 /GM 3 )^ 

is proportional to the ground meta-state of the relative motion in the hypothesis Ac R. 

If S(X, Y) = ip{[X + Y] /2)<f)(X — Y), the equation for 0, due to the spherical symme- 
try, reduces to the radial equation for x( r ) = <fr(\X — Y\). This equation has been solved 
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numerically by the algorithm obtained from the space discretization over 10 4 points of the 
equation 



1 + -iHdt/h 



u(r, t + dt) — 



1 - -iHdt/h 



u(r,t) + o(dt 2 ) 



(7) 



in the interval (— R/2, R/2), where H denotes the radial Hamiltonian for u(r,t) = 
rx(r,t). Such a procedure assures the stability of the state- vector norm during the time 
evolution and is second-order accurate [24]. 

A uniform ball of mass M = 0.38 x 10 12 m p and radius R = 4.8 x 10~ 5 cm was considered, 
with initial conditions like in Eq.(6), but for A G ~ 1.6 x lCT 6 cm replaced by A = 5.6A G . 
The solution of Eq. (3) is then obtained as the product of and the analytical solution 

\X + Y\ 2 /4 



ip([X + Y] /2,t) ocexp 



A 2 /2 + iht/M 



(8) 



of the center of meta-mass equation. 

The physical state p is evaluated by tracing out the hidden body 



p(X;X') = J dYE(X,Y)E*(X',Y). 



(9) 



After an evolution time t — 10s (~ r g ) the function g(X 1: X[) = p(X 1: 0, 0; X[, 0, 0) can 
be represented as in Fig. (4). To compare with the free evolution, in which the dynamics 
gives the usual spreading, we have also shown the corresponding function in Fig. (3). 

The final result can even be fitted by the product of two Gaussian functions: 



QiX^Xl) = exp[-(X 1 +X;) 2 /A 2 + ]exp[-(X 1 -Xi) 2 /A 2 _] 



(10) 



where A+ = 0.27R ~ 1.3 x lO^cm and A_ = 1.8 x 10 2 i? ~ 8.1 x lO^cm, while the 
free evolution, ignoring the gravitational (self-)interaction, would give the same product 
structure with A + = A_ ~ 1.3 x 10~ 5 cm (see Figs. 5 and 6). In spite of the fact that 
the fit has been performed by simple Gaussian functions, with the height and the size as 
independent parameters, the result turns out to be very accurate. 

It should be remarked that the exact rotational invariance of the initial state, which 
rather artificially leads to a degenerate final state, does not play any special role. In fact, 



due to the unitarity in the enlarged Hilbert space, modifying slightly the initial state would 
result in a slight modification of the final metastate, and ultimately of the physical state. 
Nevertheless the choice of a pure initial state corresponds to a precise measurement, which 
represents in our case a rather ideal situation. On the other hand the simplicity of the 
state constraint renders even conceptually more precise such a measurement with respect 
to the case of a system interacting with a complex environment, where the state of the 
environment after the measurement is in general a complex functional of the system state, 
to a large extent uncontrollable. In any case the present calculation has to be considered 
only as a first step in the analysis of an intrinsic nonunitary model. 



The model may also be relevant to the quantum foundations of thermodynamics, without 
resorting to the subjective and vaguely defined notion of coarse graining, by the identification 
of both ordinary entropy and black hole entropy with von Neumann entropy, i.e. with 
the entanglement entropy with hidden degrees of freedom [15]. In fact a theory that is 
intrinsically nonunitary, i.e. within which entropy can vary with time even in closed systems, 
allows in principle a much more satisfactory derivation of the thermodynamical tendency to 
equilibrium, assuming with Boltzmann an appropriate choice of the initial conditions. 

To estimate the final entropy production, consider that 8(X 2 )8(X3)g(Xi, X'- i )8{X! 2 )8{X'^) 
represents the state resulting from a measurement on p, giving X 2 = X 3 = 0, by which 
one means that the uncertainties are smaller than the typical scale of variation of p. The 
entropy of g then gives a lower bound for the entropy of p. To estimate the entropy S [g\ of 
g, we evaluate its purity 



which, for g replaced by its analytical fit (Eq. 10), gives Trg 2 = A_/A + ~ 6. x 10 2 . If, for 
simplicity, we consider the corresponding ensemble as one of N equiprobable states, then 



IV. ENTROPY ESTIMATION AND FINAL REMARKS. 




(11) 
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Tr ^ = Ei^ = ^ iV ~ 17 ' S0~/r B logl7. (12) 

3=1 

If we approximate p as the direct product of three equivalent g, namely we omit the entan- 
glement between the three Cartesian coordinates, which is absent in the initial pure state 
and would stay so if the potential in Eq. (4) were replaced by its quadratic approximation, 
then the total entropy is S [p] = 3S[g\, corresponding to N 3 equiprobable states. 

Independently of any approximation p(X; X'), as the kernel of a compact positive semi- 
definite Hermitian operator of unit trace [25], can be diagonalized as 

p(X;X') = Y^MAXWAX'); p 3 > 0; J2P3 = 1; (^#fc> = 8jk, (13) 

3 j 

where the above approximate estimate makes us expect that — pj logpj ~ 3 log 17. This 
result has to be compared with a value ~ 31og[A + /Ac] ~ 3 log 8 corresponding to a naive 
counting where the orthogonal states ipj above are assumed localized and approximately 
non overlapping. This small discrepancy is not surprising, as these states are expected to 
include contributions from several low lying bound meta-states, so that they are orthogonal 
in spite of the overlapping of their probability densities, due to their space oscillations. 

The qualitative agreement between our two estimates of entropy, the one corresponding 
to the computed density matrix and the other corresponding to the approximation of the 
mixed state by means of an ensemble of equiprobable localized states which occupy the 
volume roughly occupied by the density p(X; X), makes it natural to assume that the states 
tpj diagonalizing p are localized. To be more precise, this would be strictly true, without 
ambiguities, after breaking the exact rotational invariance, which, as mentioned, is expected 
to introduce artificial degeneracies in the density operator. 

The present result strengthens our confidence in the most relevant peculiarities of the 
localization phenomenology ensuing from the model, which make it in principle distinguish- 
able both from the other proposed models [7,8] and possibly from the competing action of 
the environment-induced decoherence [22]. In particular the model presents a sharp thresh- 
old, below which localization is practically absent, and a localization time r g oc M -5 / 3 
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rapidly decreasing, as the mass is increased. Furthermore the threshold mass M t oc p 1 / 10 is 
remarkably robust with respect to mass density variations. 
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Captions : 

Fig. 1. Initial radial wave function x( r ')> with r' = r/R, for the internal state of the 
two meta-bodies. 

Fig. 2. Radial wave function x( r ')> with r' = r/R, for the internal state of the two 
meta-bodies at the final simulation time t — 10 s. 

Fig. 3. The modulus of the reduced density matrix \p(X,X')\ in the absence of gravity 
is shown in the plane Xi,X[ with X 2 = X' 2 = X 3 = X^ = 0, at the evolution time t = 10 s. 

Fig. 4. Modulus of the reduced density matrix \p(X,X')\ with the inclusion of gravity 
in the plane X±,X[ with X 2 = X' 2 = X 3 = X% = at the evolution time t — 10 s. 

Fig. 5. A longitudinal Gaussian profile corresponding to the function \p(X, X) \ has been 
superimposed to the array of points obtained from the numerical simulation. Lengths X on 
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the x-axis are measured in units of the radius R. 

Fig. 6. Gaussian transverse cross-sections corresponding to the function 
\p(X, —X + k * 0.08/2) | , k — 0, 1, 2, 3, 4 have been superimposed to the arrays of points 
obtained from the numerical simulation. The dashed line represents the longitudinal cross 
section. Lengths X on the x-axis are measured in units of 10~ 2 R. 
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